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CjT)! Abstract 



In analogy with classical electromagnetic theory, where one determines the total 
charge and both electric and magnetic multipole moments of a source from certain 
surface integrals of the asymptotic (or far) fields, it has been known for many years - 
from the work of Hermann Bondi - that energy and momentum of gravitational sources 
could be determined by similar integrals of the asymptotic Weyl tensor. Recently we 
observed that there were certain overlooked structures, defined at future null infinity, 
that allowed one to determine (or define) further properties of both electromagnetic 
and gravitating sources. These structures, families of complex 'slices' or 'cuts' of 
Penrose's 3 + , are referred to as Universal Cut Functions, (UCF). In particular, one 
can define from these structures a (complex) center of mass (and center of charge) 
and its equations of motion - with rather surprising consequences. It appears as if 
these asymptotic structures contain in their imaginary part, a well defined total spin- 
angular momentum of the source. We apply these ideas to the type II algebraically 
special metrics, both twisting and twist-free. 
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I. Introduction 



In the present work we return to the study of an old subject p 19 | 3 | /t | namely prop- 
erties of asymptotically flat Einstein space-times, both vacuum and with a Maxwell 
field. More specifically we explore some previously overlooked properties of such 
space-times. The point of view is however rather new. 

The basic idea is that for each asymptotically flat vacuum solution there is a 
certain 'structure', a unique function on null infinity, 3 + , that has been known in 
various contexts, but whose several unusual properties have largely been overlooked. 
This function, which we shall refer to as the Universal Cut Function (UCF), carries 
virtually all the asymptotic information, as well as much of the physical content of 
its associated solution. In addition, it unifies several apparently disparate ideas. In 
the case of the Einstein-Maxwell fields, in general there are two such functions, one 
for GR and one for the Maxwell field - however we will consider only the degenerate 
case where the two function are the same. 

We begin with the conventional view of 3 + ; it is a three-dimensional manifold, 
S 2 xR, coordinatized by the standard Bondi coordinates (u r ,(, C), with complex stere- 
ographic coordinates ((, () on the S 2 and u r on the R part. (The subscript V denotes 
the retarded time.) We part from the conventional point of view by allowing the com- 
plexification of 3 + , i.e., we consider 3^ with u T now allowed to take complex values 
close to the real and (£, () going to (£, () with ( close to (. It should be emphasized 
that this complexification of 3 + is a formal device that has many interesting conse- 
quences. Nevertheless we are dealing with real space-times and real solutions to the 
Einstein equations. 

A significant notational problem should be addressed now. The time variable u r 
used above is, in flat space the retarded time t — r. Many years ago,^ in order to keep 
the symmetry of many of the equations in the spin-coefficient formalism this u r was 
replaced by u, via 

u r = V2u. (1) 

Later to see the physical consequences of our results we restore the retarded time u r 
and, in addition, explicitly display the velocity of light, c. 

In Sec. II we review some background material. Much of it, though in the liter- 
ature, is still quite unfamiliar to most readers. To help overcome this unfamiliarity, 
we present in Sec. Ill a detailed discussion, in the simplified case of flat space-times, 
of the ideas to be used later in the asymptotically flat cases. Sec. IV is devoted 
to a discussion of these ideas for general asymptotically flat space-times while Sec. 
V is used to specialize the discussion to the algebraically special type II Einstein- 
Maxwell metrics. In Sec. V and the remaining sections, we first see how the new 
ideas interact with the classical Bondi energy-momentum four-vector, giving it new 
properties and even fresh predictions. Then we see how it interacts with the Bondi 
energy-momentum loss equation, giving it surprising new interpretations that include 
equations of motion for the center of mass motion, even containing radiation reaction 
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forces. It even strongly suggests a general relativistic definition of spin-angular mo- 
mentum. Finally in Sec. VIII, the conclusion, we discuss how in principle one could 
try to observe the predictions described earlier. 

II. Background 

The UCF, which had been introduced earlier in different contexts] \4MU j s described 
as follows: It is to be a one-complex parameter family of 'cuts' - complex 2-surfaces 
- of 3q written as 

u = X(t,C,Q, (2) 

with (in general) a complex r parametrizing the family. [In important special cases, 
we could have a real one parameter family of real cuts of 3 + , so that u = X(t,(,C) was 
a real function of a real r. In fact, in this case, Eq.(T5]) could be treated as a coordinate 
transformation from u to r. The real coordinates (r, £, Q have been referred to NU 
coordinates.] We also assume that Eq.j2j) can be inverted, i.e., solved for 

r = TO, CO- (3) 

For real u and C = C this obviously defines a complex r such that the u in Eq.((2]) is 
real. {Aside: Later we will see that T for real u and ( = ( is a CR function for a CR 
structure obtained from X(r, (,()■} 

From X(r, (,() we define two functions, V and L, (that have been extensively 
used in GR in different contexts, with different coordinate systems,^ by 

V = d T X = X'(T,CO, (4) 
L = 5 (T) X(r,C,C), (5) 

with 5( T ) meaning the edth operator with r held constant. Here, these functions 
have a common origin in the UCF. An important point is that, in both V and L, 
the r could be eliminated via Eq.((3|) and they could be treated as functions on the 
real 3 + . Both have simple geometric meanings: V = X'(t,(,() is the V separation 
between neighboring constant r cuts while L(u, (, () is a stereographic 'angle' field 
on real 3 + , on the sphere of past null directions, that determines the directions of an 
asymptotically shear-free null geodesic congruence going into the interior. 

In the case of algebraically special type II space-times, V(r, £, C) is the basic 
variable in the Robinson- Trautman metrics^ with real r, while L(u,(,C) is the basic 
variable for the twisting type II metrics. 

• Note that we use (') = d T and (') = d u 

Of primary importance to us is the fact that X and hence L and V are uniquely 
determined (up to gauge transformations - BMS transformations) in any given (non- 
flat) asymptotically flat space-time. They are determined by the appropriate choice of 
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a shear-free (for fiat or algebraically special space-times) or an asymptotically shear- 
free (for asymptotically flat space-time) null geodesic congruence. 

For flat space-time, using the variables (u, CO? X and L are determined by first 
solving the equation for shear-free null geodesic congruences,^ 

<5 {u) L{u,(,0 + LL- = 0, (6) 
for L(u, C,C) followed by solving the so-called CR equation 

5 (u) T + L(n,C,C)T=0, (7) 
to determine r = T(w, (, (). Finally by inverting T(u, (, (), we obtain u = X(t, (, (). 



• 5(„) is the eth operator with u held constant. 

An alternate and easier method to find X and L , is to immediately use the 
variables (r, C, C)> an d solve the equation 

af T) x(r,c,C) = o, (8) 

for u = X(r, C, C) an d then find L(r, (, Q by 

L(t,C,C) = 3 (t) X(t,C,C), 

noticing that 

5 (T) L(r,C,C)=0. 

These results follow by implicit differentiation of Eq.(T51) using Eq. ([6]) . In the flat 
space case there is considerable freedom in the choice of solution for X and L that is 
considerably narrowed by the requirement of regularity of the function L on 5* 2 . The 
geometric meaning of 'regularity' is that all the null geodesies are transverse to 3 + , 
i.e. none of the geodesies lie on 3 + . The final freedom, which follows from Eq.(j8]), is 
that of a free complex analytic word-line in complex Minkowski space. We return in 
more detail to this issue in Sec. II. 

For asymptotically flat space-times with a given Bondi asymptotic shear, a = 
a(u, C, (), generalizing the flat space treatment, the UCF is determined by first solving 

<5 {u) L{u,(,0 + LL = a{u,(,0, (9) 
for L(u, (, (), followed by solving for T the CR equation 

a H T + L( M ,C,C)T=0. (10) 

By inverting the solution, Eq.([3|), we obtain X(r, £, C)- Again there is considerable 
freedom in the solution which is first narrowed again using regularity on S 2 , then 
made unique by certain algebraic conditions on the Weyl tensor components. 1 11 1 12 1 1 -^ 
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An alternative method of finding the UCF from a given a(u, (, () is to solve the 
PDE, the so-called good-cut equation] 14 1 15 1 



5 2 z(c,C) = ^,C,C). (ii) 

The solutions, which are known^ED i depend on four complex parameters, say z a } 
can be written as 

Z = Z(z a X,0- (12) 

When z a are taken as arbitrary complex analytic functions of r, i.e., a complex world- 
line in the parameter space, 

z a = C(r), (13) 

we have 

u = X(r,C,0 = Z(C(r),C,0, (14) 

which is almost, but not quite, the UCF. 
At this point we have: 

Theorem 1 Every transverse (regular) shear free or asymptotically shear-free null 
geodesic congruence in an asymptotically flat space-time is generated by a complex 
world-line in H-space. 

Finally an algebraic condition on the Weyl tensor components determine d 11 1 12 ' 13 ^ 
a unique complex world-line and thus yields the UCF. 

One knows from electrodynamics that many properties of the sources (charge and 
current distributions) of an asymptotically vanishing electromagnetic radiation field 
can be determined or at least studied from the asymptotic field itself, e.g. the total 
charge and both the electric and magnetic multipole moments. In GR one tries to do a 
similar thing. For example we have the Bondi energy-momentum defined from certain 
asymptotic Weyl tensor components and many attempts to define multipole moments 
and angular momentum from other Weyl tensor components and the asymptotic 
shear. 

Our contention is that, aside from the energy-momentum four-vector, many or 
most of the asymptotically obtained physical quantities are hidden in the UCF. And 
even the properties and description of the energy-momentum four-vector is greatly 
clarified by the UCF. In order to extract this information from the UCF , using its 
assumed regularity, we expand it in spherical harmonics, writing it as 

U = X(T,(,0 = ^ lm (T)Y lm (CX), 

or, using the spherical tensor harmonics^ as, 

u = X(t,C,C) = ^e°(r)F (C,C) - \?(t)Y°(M)+P(t)Y&{M) + (15) 



5 



where the components ^ lm (r) or the (£°(t), £*(t), £^(t), ) are ) m general, complex 

functions of a complex time r, and are to be given physical meaning. Eventually, the r 
must be replaced, via Eq.((3|), by the real u. The four components of £ lm )-v^(£°, C) 
are to be identified with a complex position vector, the real part with the center of 
mass, and the imaginary part with spin-angular momentum/unit mass. The complex 
£ 2m -vv- (? 3 are associated with mass and spin-quadrupole moments while the higher 
harmonic terms are related to the higher moments, the real parts to the mass moments 
and the imaginary parts to the spin-moments. 

In addition to its role as describing physical quantities, the UCF has a pretty 
geometric/group theoretical meaning. 

One knows from Bondi's and Penrose's early work concerning the asymptotic 
symmetrie s ' 18 1 19 1 7 1 8 - 1 that the freedom in the choice of the Bondi coordinate system is 
what is known as the BMS group. In addition to the rotation subgroup it contains 
the supertranslation subgroup acting on 3 + and given by 

u = u' + a(C,(), 

which can be extended into the complex by 

u = u' + a(CX). 

The function £), which can be expanded in spherical harmonics, contains the 
(complex) Poincare translation subgroup defined by the expansion in only the I = 
and I = 1 harmonics, 

«(C,0 = M(C,C)+^(C,0- (16) 

If we choose the starting value u' = and consider a one parameter family of su- 
pertranslations, a(C, C) =^ a ( T XX)i we can consider the UCF as a one complex pa- 
rameter family of supertranslations and the I = and I = 1 part as a one-parameter 
family of complex Poincare transformations. 

In other words, the complex four-vector world-line £ a (r) = (£°(r), ^( T )) can be 
interpreted as a world-line in the (complex) Poincare translation subgroup. We can 
chose t so that the velocity vector of the world-line has a Lorentzian norm of one. 
{An alternative interpretation is to take the world-line in H-space and use the H-space 
norm.} 

The remainder of this work is organized as follows: In Sec. Ill, the application 
of the UCF to Minkowski space is described and the two different flat space forms of 
the metric, those associated with L(u, C, C) an d with V(r, £, C), a re worked out. In 
Sec. IV we give a brief outline how the UCF is used in general asymptotically flat 
space-times while in Sees. V-VII we apply these idea to the algebraically special type 
II metrics. Sec. VIII is devoted to a discussion. 
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III. The UCF and Flat Space 



Since many of the ideas developed here are unfamiliar to large portions of the 
relativity community we believed it would be of value to first describe them in the 
simpler context of flat space-time where the functions L and V have made frequent 



earlier appearances^^ and have relatively simple meaning. 

The basic idea here is to first describe arbitrary shear-free null geodesic congru- 
ences in flat space and then to specialize them to the case of regular congruences, 
i.e., congruences that are completely transverse to 3 + \ In general these congruences 
have twist and the basic variable of choice is the L. Twist-free congruences are more 
easily treated with the variable V. The main distinction between the flat space and 
the asymptotically flat space properties of L and V is that in the spherical harmonic 
expansions only the lowest harmonics, i.e., / = and 1, appear in the flat space while 
all harmonics can appear for the curved space. 

We begin by giving an analytic description of any arbitrary null geodesic con- 
gruence, (with l a as tangent vector), in Minkowski space^ . Using x a as standard 
Minkowski-space coordinates, L(u, C, £) an arbitrary complex function of the param- 
eters (u, C, () which label individual members of the null geodesic congruence and r 
as the affine parameter along each null geodesic, we find that 

x a = Ur t a - Lm a - Lm a + (r - r )l a , (17) 

x a = u r K ' - Lfn a - Lm a + (r - r )l a , (18) 

a/2 



u r 

with the null tetrad 



x a = u{l a + n a ) - Lm a - Lm a + (r - r )/ a , (19) 
: V2u, (20) 



r = + CC, C + C, -<(C-C),-i + C0; (2i) 



2P 



V2, 



m a = 3P = J^(0,l-C 2 ,-z(l + f), 2C), 

rff = sr = ^|(o,i-c 2 , *(i + c 2 ),2C), 
n a = -C + C, i(C- 0,1 -co, 

la — I — 

e - (22) 

P = 1 + CC (23) 

For the simplification of the complex divergence of the congruence the function 
r is most often chosen as 

r = --{tiL + LL +KL + LL-}. (24) 
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Our main interest lies in the class of regular null geodesic congruences with a 
vanishing shear. This condition is achieved by first imposing the shear-free condition 
on L, i.e., the differential condition^ 

5L + LL = 0, (25) 

and then requiring solutions to be regular functions (expandable in spherical harmon- 
ics) on the (C, C) sphere. To solve this equation we introduce the (complex) auxiliary 
function, {a CR function} 

t = T(u,(,(), (26) 

by 

ar + LT=o. (27) 

By using the inverse function to r = T(u, (, (), i.e., the UCF 

u = X(t,(X), (28) 
and implicit differentiation, L can be written as 

L = a (T) X(r,C,C), (29) 

and Eq. ([25]) becomes 

5? T) X = 0. (30) 
The general regular solution, L(u, £, £), to Eq.(j25l) is then given in the parametric 

form 

L = e(r)m a {CC), (31) 

u = e(r)l a (C,0, (32) 

where ^(t) are four arbitrary complex functions of r. Since Eq.jH]) is invariant under 
T* = F(T), we can take r so that 

v a {r) = e(ry = o T e, 

is a unit Lorentzian four- vector, v a v a = 1. Using the spin-s tensor notationpEl 

i i 

a ~ l ^2' 2 uh 

m a = (0,^), 

Eqs.([3l]) and ([32]) become 

L = C(t)Y', (33) 

-0/ 



u = ^-\C(r)Yl (34) 
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Remark 1 We comment again that L is a complex function that describes a stereo- 
graphic 'angle' field on real 3 + , i.e., a direction field that determines a null geodesic 
congruence. It is generated by the UCF that is determined by an arbitrary complex 
world-line in complex Minkowski space 

z a = C(r). 

When L from Eqs. and (34\ ) is used in Eq. ( TTPp . we obtain the explicit form for 
any regular shear-free null geodesic congruence. 

Later, in more general situations, we will use this form, Eqs.(l33l and (1341) , for 
the first two harmonics. Curvatures come from the higher harmonics. 

An alternative meaning of Eq. (fT9l) is that it defines a coordinate transformation 
from the natural coordinates of the congruence (u, r, (, Q to the Minkowski coordi- 
nates, x a . One can in a straightforward way find the flat metric in these coordinates. 

From Eq. (fT9l . calculating the differentials we find 

d x a = (1 - r )l a du + n a du + l a dr - L'm a du - Lm a du (35) 
+P- I [(r - r - 5L)m a - ~oTm a - Ln a + (L - ~dr )l a }d( 
+P- I [(r - r - 5L)m a - 5Lm a - Ln a + (L — Sr )l a }d(. 

Since the set (/ = l a dx a ,n = n a dx a ,m = m a dx a ,m = m a dx a ) define a null one-form 
system the metric can be written as 

V = Vabdx a dx b = 2{ln — mm). (36) 

Any other set of one-forms obtained from the first set by a Lorentz transformation 
yields the same metric. Choosing a new null tetrad via the null rotation 

I* = I, 

m* = m — L'l, 

m* = m — L I, 

n* = n — L m — L'm + L L'l, 

we obtain, by a direct (lengthy) calculation, using Eq.(l25il. the simplest form of the 
tetrad, namely 

I* = dx a l a = du- P^iLdC + Ld(), (37) 
m* = -P- 1 [r-iE]d( = P- I p- I d(, 
m* = -P- 1 [r + ii:}d( = P^p^dC, 

n* = dr + p- 1 {K + Lp~ 1 )dC + (K + L p- l )d~Q + (1 + i{5L + "5L' + LL" +LP })1, 
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where 

P = 7^F> (38) 

zE = -r -3L -LL = -{$L + LL -bL -LL}, (39) 

K = L + Lar + ^{5 2 I + L 2 I" + L[5l-5L + Lr-lL]}, (40) 

with p being the complex divergence of the congruence and E its twist. The metric 
is then given by. 

7] = 2(/*n* -m*rff). (41) 

This form of the metric is very close to that the Type II twisting metricsP 
If the world-line £ a (r) was described by a real function of a real r, though Eqs.(l37l) 
and (14B could still be used {with now E = 0}, there however is a simpler set of 
coordinates, namely (r, r, CjC)j associated with the non-twisting congruence, that 
can be used.^ They are given by 

x a = C(r) + ^l a (U), (42) 

where 

u = C(r)l a = X(rXX), (43) 

v = r = a = A = |-^. 

By a computation similar to that with the complex curve we obtain the flat metric 
in (r, r, (, () coordinates, 

ds 2 = r] ab dx a dx b = (1 - 2^r)dr 2 + 2drdr - r 2 ^|, (44) 

which is very similar in form to the type II (non-twisting) Robinson-Trautman met- 
rics.™ 

As we pointed out, in flat space the UCF is not unique. It is determined, in 
general, by an arbitrary complex world-line in complex Minkowski space. However in 
the presence of an asymptotically vanishing Maxwell field (with non-vanishing total 
charge) one can construct, from the Maxwell field itself, a unique complex world- 
line (the complex center of charge)^ and from it a unique UCF can be found. For 
example, the Lienard-Wiechert charged world-line defines^^the real £ a ( r )- 

In the case of asymptotically flat vacuum space-times, properties of the Weyl 
tensor are used to define the world-lin e 1 11 1 12 1 13 1 1 ^ and the unique UCF. In the case of 
asymptotically flat Einstein-Maxwell space-times there are in general two different 
complex world-lines, determined respectively by the properties of the Weyl tensor 
and the Maxwell field. As we said earlier, we here consider only the case where the 
two lines coincide. 
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IV. Asymptotically Flat Space-Times 



In this section we outline - leaving out many of the details - how the UCF is 
determined in asymptotically flat space-times. 

For asymptotically flat space-times, the basic radiation data is the Bondi asymp- 
totic shear, o~(u, C, C)i a complex function given on 3 + which, as we mentioned earlier, 
we assume can be analytically extended a short way into Jq, and hence written as 

The a(u, (, () describes the shear of the null geodesic congruence that are the 
null generators of the Bondi null surfaces given by u = constant in the neighborhood 
of 3 + . Any new asymptotic null geodesic congruence can be obtained by a Lorentz 
transformation (specifically a null rotation around the null generators of 3q) at J + . 
If l a are the tangent vectors to the Bondi congruence the new tangent vectors, l* a , in 
the neighborhood of 3 + are given 

l* a = l a B -^m a B + ^-m a B + 0{r- 2 ), (45) 

m* a = m * B --n a B + 0{r- 2 ), 
n* a = n%, 

with L(u,(X), a t this point, an arbitrary complex function on 3 + . The asymptotic 
shear of the new congruence is given by™ 

a* = o — 5( U )L — LL' . (46) 

By requiring that the new congruence have vanishing asymptotically shear, we see 
that L(u,(,C) must satisfy the differential equation, Eq.((9]), 

Z (v) L(u,<Z,C)+LL- = a(u,C,0, (47) 

where now we are considering the extension into the complex. The technique for 
studying properties of solutions of Eq. (11711 is quite similar to what was done in the 
Minkowski space case, Eq. (l25l) . 

Again we introduce the new complex auxiliary function 

r = T(u,(,Q, (48) 

by 

r- 9 H T (49) 



T- 7 
and its inverse 

u = X(t,(,(), (50) 



11 



which eventually leads to the UCF. {Note that we have the useful gauge freedom 
r* = T*(u,(,() = F(T(u,(X)) in the determination of T(u, £,£)■} By implicit dif- 
ferentiation of Eq. (l50l) . i.e., 

1 = X'T\ 

= a (T) X + X'5 (u) T, 

we find the implicit form of L(u, C, C) 

L = a (T) X(r,C,C), (51) 
u = X(t,C,0- (52) 

Continuing with the change of variable, Eq.(l47l) becomes 

b 2 (T) X = a(X,CX), (53) 

a well studied equation, frequently referred to as the 'good-cut' equation. 1 14 1 15 1 Its 
solution, which depend on four complex parameters, z a (defining H-space) is written 

as 

x = x(z a ,(,0- 

Now by taking an arbitrary world-line in the parameter space, z a = £ a (r), we have 
the parametric form of the solution to Eq.(T47l ). 

u = X(r,C,C)=X(e a (r),C,C), 
L = 5 (T) X(r,C,C). 

It is important to notice {see Eq.fl3"4l)} that in the harmonic expansion of X(r,(,() 
we can take the first four harmonic coefficients, / = (0, 1), as the £ a (r), i.e., write 

u = C( T )L + harmonics (I > 2), (54) 
to 1 

u = -^=(t) - -C(r)Y^ + harmonics {I > 2), 
v 2 2 

and that 

v° 1 

V = X' = — =(r) - -v*(r)yS + harmonicsil > 2). (55) 

We see that in both the flat-space case and for asymptotically flat space-times 
the freedom in choosing a regular shear-free or asymptotically shear- free null geodesic 
congruence lies in the arbitrary choice of a complex world-line in the parameter space, 
i.e., complex Minkowski space or H-space - or via the alternative interpretation, the 
space of complex Poincare translations. 

The issue is: can one uniquely determine this world-line'! 
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In Minkowski space the answer is no, except when there is an asymptotically 
vanishing charged Maxwell field. Using the null rotation, Eq.fl45l). and requiring the 
vanishing of the coefficient of the / = 1 harmonic of the Maxwell tetrad field, 4>l°, 
(the complex dipole term), one determines a unique world-line, the 'complex center 
of charge' pi 

For asymptotically flat spaces, with a similar null rotation, and then requiring the 
vanishing of the I = 1 harmonic of the Weyl tensor component, ty®, i.e., the complex 
mass dipole terms, the world-line 

z a = C(r), 

is again determined) 11 * 12 ^ the details for the determination of the unique complex 
curve have been described earlier^^l and a further discussion would take us far afield, 
we will not pursue this issue here. 

Remark 2 We note that the information in the Bondi shear begins with the I = 
2, i.e., quadrupole, harmonic while the UCF, in addition to its higher harmonics, 
I ^ 2, also contains the information in the I = 0, 1 harmonics obtained from the 
unique world-line. In other words while the UCF uniquely determines the Bondi 
shear, Eg. (53) . the Bondi shear only determines the higher harmonics of the UCF. 
In essence the new structures we have been describing here are just these I = 0, 1 
harmonics of the UCF or equivalently, the complex world-line. Later we will see the 
physical content contained in this new information. 

As an aside we note that in the very special case, when the unique asymptotically 
shear free congruence is also twist-free, the UCF, u = X(r, C, (), will be a real function 
of a real r. It could then be interpreted as simply the coordinate transformation 
between Bondi time and NU time. 



Remark 3 To end this section, we describe the natural relationship our UCF has 
with the elegant mathematical idea of a CR structur ^ 2 ^ 25 ^ 27 ^ . The UCF defines 
or induces a CR structure on the real J + in the following fashion: from the UCF 
we obtain L(u,C,C)? as above, but now restrict (u, C, C) to the real 3 + . Then define 
the CR structure by (the equivalence class of) the three one-forms, one real and the 
complex conjugate pair: 

1* = du ^— dC — dC, (56) 

1 + CC 1 + CC 

* dC _* dC 

m = -=■, m 



i + CC i + CC 

Writing ( = x + iy, the functions r = T{u, x + iy,x — iy) and ( = x — iy are 
the CR functions defining the embedding of 3 + into C 2 = {t,(}, both satisfying the 
CR equation, i.e., d( u )T + LT' = 0, etc., which was originally our defining equation 
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for T. It is natural since the complex conjugate pair ( m*, m* ) are the dual covectors 
to the u = constant surfaces of 3 + and I* is the real one-form determining the null 
direction of the new asymptotically shear-free null geodesic congruence. 

Though we will not make further direct use of the CR structure we note that it 
is the CR structure that lies at the basis of much of the present work. The function 
r = T(u, C, C) essentially defines the CR structure and its inverse, u = X(t,(,(), 
yields the physical content of the work. 

V. Type II Metrics and the UCF 

A. Non-Twisting:The Robinson- Trautman Case 

Though our issue here is the analysis of the type II twisting metrics - with and 
without a Maxwell field - we nevertheless want to begin with a brief discussion of the 
non-twisting case, i.e., the Robinson-Trautman metrics. 

Since Robinson-Trautman (RT) metrics have been extensively studied and writ- 
ten about EHHQ] we w jjj con fi ne ourselves to the discussion of their relationship with 
the UCF. 

Rather than beginning with the differential equations for the determination of the 
Robinson-Trautman (RT) metrics, we reverse the proceedings. We will assume that 
we have the general solution to the RT equations and that furthermore we explicitly 
know the unique UCF, 

u = X rt (t,(X), (57) 

which in this case is real with real r. The (r, (, () are the RT/NU coordinates on J + , 
with the UCF, Eq. (l57l) . having the added meaning of being the coordinate transfor- 
mation to Bondi coordinates. The function 

V rt = X' rt (t,U), 
then is the solution to the RT equations 

5 (T) x = 0, (58) 

X ' ~ - V^&I^Vrt - V^b 2 {T) V RT ■ b 2 {T) V RT } = 0, (59) 

Vrt 

X = -V^f. (60) 

where ^ is a Weyl tensor component. From Eq. (l58|) we see that x depends only on 
r. Using the freedom of time change, r = /(r*), x could be made constant. However 
this freedom was already used to make v a {r) = £ a/ into a unit vector so that x 
remains a function of r. The time evolution, Eql59l. is in fact equivalent to the Bondi 
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energy-momentum loss equation. When Eq.fj59l) is transformed to Bondi coordinates 
via Eq. (l57l) it becomes 

*• = (j(wy, 

with the mass aspect given by 

= $f + 2L6(a)- + L 2 (a)~ + 5 2 a + a(a) = (61) 

The L and a are directly calculated from the UCF as described earlier. 

Though it is virtually impossible to exactly perform the coordinate transforma- 
tion to Bondi coordinates, one can do so within an approximation scheme. Using a 
second order expansion, essentially in v/c, one finds from the mass aspect that the 
Bondi mass and momentum can be partially expressed in terms of the world-line and 
take on the surprising familiar kinematically form 

v i2 v i/£i 

M = M [1 + — — H f-1 + quadrupole terms, (62) 

2c 2 & 

P l = M v l + quadrupole terms, 
where we see the approximation to the relativistic mass 



with the strange term Mq—§- appearing to be the "work done" moving the "particle" 
to the "position . 

The Bondi mass-momentum loss equations, to second order, become the equa- 
tions of motion for the world-line. 

W = --^GQ ij '"Q ijm , (63) 
be' 

Mv" = v l -^jQ jk '"Q jk '", (64) 

DC 



where Q is defined from the I = 2 coefficient in the harmonic expansion of X, Eq.ffT5il 
by 

( " = ^ Q '"' (65 » 

These results, when generalized to the Robinson- Trautman-Maxwell equations,^ 
contain further surprising results, e.g., Eq.(l64l) contains, among other terms, the 
classical radiation reaction force given by 



2q 2 
3c 



See also Eq. (11381) below. 
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B. Type II Twisting Metrics 



1. Solving the field equations 

Our primary objective in this work is to analyze and look for the physical meaning 
of the different geometric variables that arise in the vacuum and Einstein-Maxwell 
type II twisting metrics. In particular we want to see how the UCF enters into 
the discussion. The relevant final differential equations to determine the metrics have 
been known and studied for many years J EEEEEQl Though a great deal of mathematical 
insight, specifically the discovery and understanding of the associated CR structure, 
was developed, aside from the Kerr metric and the charged Kerr metric, general 
physical understanding of these metrics seems to be missing. Here, with the help of 
the UCF, this lack of understanding will be addressed. 

Acknowledging the non-linearities and the difficulties in finding exact results, 
we approach the problem of physical interpretation using, in addition to the UCF, 
the two approximations: (1) In the spherical harmonic expansions only terms with 
/ = (0,1,2) harmonics will be retained and (2) calculations will be done to second 
order in v/c. 

The Type II coupled differential equations for the independent variables (L, $1°, 4$ , 
in terms of KCCP 1 ™ 51 

801° + 2L'<t>f + L{4>fy = 0, (66) 
(0f)+80;°+(L0; )- = 0, (67) 



= 8^ O + L(^°)' + 3LV;°-2A;0*V2 , (68) 

vD- = a yy + k<p* 2 %°, (69) 

# = ^ = ^° + 2Ld(a)- + L 2 (cf)-- + d 2 cr + a(cf)-, (70) 

with 

cr = 8L + LL . (71) 

Eq. (i70l defines the mass aspect with the requirement that it be real while Eq.(l7Til 
defines the Bondi shear a from the L. Eq.(l69l is the evolution equation for the 
Bondi mass aspect with the / = (0, 1) terms determining the energy and momentum 
loss. Basically these equations have the following structure: from the two Maxwell 
equations, Eqs.(l66l) and ( l67l) , both (f>*® and 02° will be determined in terms of the L 
and a constant of integration, q, the charge. Eq. (l68l determines ip*® (or \I>) in terms 
of L and a single function T(u). The function L(u,(X) could then be considered as 
our basic variable which would be determined by Eq. (l69l) . We however will consider 
something different. Our point of view is that L is a derived function, coming from 
the UCF as was described earlier. Since it is the information encoded in the UCF 
that we want to use, we must reexpress the differential equations in terms of the UCF. 
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Our first task will be to transform three of the field equations, i.e., Eqs.(l66l). (l67j) and 
681 ). to a new very useful form. 

Assuming that we know the UCF, u = X(t,(,()i then and its inverse, r = 

T(u, (, () and T' = d u T contain the information from X. We define T(u, (, (), C> C)> $2(u, C C) 
by 

r 2 ° = (TfT, (72) 

#° = (T) 2 $ 1; (73) 

<t>f = T'$ 2 . (74) 



When these expressions are substituted into Eqs.(l66i). ( 1671) and [68]), remembering 

that 

L = (75) 

their new form becomes 

8$i-^$i = 0, (76) 
5$ 2 -^$ 2 = -(TT^iiTf}-, (77) 

[ar-^x] = 2^$^. (78) 

It is these three equations that we will solve within our approximation scheme. 
Remark 4 Eg. Fffl) can be integrated immediately by 

$! = $!(:?>, c,0), (79) 

with $1 an arbitrary function of T. This follows from regularity on the sphere 
and the fact that Eg. f/h)) defines a CR function. 

The main difficulty lies in finding or calculating the T(u, (, £) from the assumed 
form of the UCF, X(t,(,(), i.e., from Eq.(l54j), 

u = X(r,C,C) = ^°0W(C,C) - \e(r)Y^CX)+e j (r)Y^(C,~0 + •-, (80) 
All the relevant functions are assumed to be analytically extended a short way into 

T(u, C, C) can be approximated in the following manner: From the unit velocity 
vector v a {r) = £ a (r)' = (v°(t),v 1 (t)) with a slow motion approximation 

v \ T ) = ^/l + v i2 = l + \v i \ (81) 
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we have 



e°(r) = J v°(r)dr = r + 5£°(r) = r + ^jv t2 dr + & 
The Eq.(l80l) can then be written as 



or as 



2 

2 

r = ^ + F(r,C,C) ; 



^ = T+5e(T)-^e(T)Y°(cx)+V2e j (T)Y° j ((,o + ---, 

r = ^-^^r) + ^(r)F°(C ! C)-^(r)r4.(C,C) + ... 



with 



F(r, CO = -6?{t) + V^OWiCCC) " ^(^&(C,C) + • 



72 
2 

The idea is to iterate Eq. (j83|) . 
Taking 



r = \f2u = w r , 



as the zeroth iterate (see Eq.(pQ)), we have for the first iterate that 
r = V2u+F(V2u,C,0, 

r = Ur + ^CMY^CO-V^C'MY^CX), 
and for the second iterate 

T = T{u,(,0 =u r + F(u r + F(u r ,(X)XX) ~ u r + F + Fd Ur F. 
We then find that 

X" = v^T = v^(l + F' + (F') 2 + FF") , 
or expanded and using the Clebsch-Gordon products, 

- ^K'e + v 3 'C j + 2«*V}]yg 

+^K'£ fc + - QV2v ki + ^[£*V + v kl v il }}Y 2 ° kl ). 
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Remark 5 We have used the derivative notion d u = (') and d T = ( ). To avoid a 
plethora of notations we also use d Ur = f) since the functions ofu r (after expansions) 
are always functionally the same as those ofr. Which variable, r oru r , is being used 
is clear from the context. Furthermore we have d u = (') = \/2d Ur = v2 ('). All the 
harmonic coefficients in the remainder of this section are functions only of u r . 

Functions just of T, as for example the $1 of Eq.(l79il. can be easily approximated 

as: 

$(T) = $(u r + F(u r ,(,() + Fd Ur F) (90) 
= $(u r ) + (F + FF')$', 

or, when $ is already first order, 

$(T) = $(u r ) +F&. (91) 

We will take $i(T), with q = constant, to have the form 

$ 1 = $ 1 (T) = | + 5$ 1 (T). (92) 

The quantities T, $x and $2 are now assumed to have the harmonic expansion 
given by 

$! = <s> 1 + <s>\Y° + ^Y° j , (93) 

$ 2 = $>\Y u l + c&f^T 1 , (94) 
T = T° + T*y x ° + T i3 V 2 ° 3 -. (95) 

By substituting these into the field equations, Eqs.jjH), ((771) and ([78]), (or Eq.([9l]) 
for $1) and using Eqs. (l87l) and (I89al) . after very lengthy calculations with frequent 
use of the Clebsch-Gordon expansions we obtain the solution to the new version of 
the Maxwell equations 

0f = (T') 2 $i, (96) 
SiCO = | + 5$ 1 (T) = |(l-^ 2 ), (97) 



and 



(T)$ 2 , (98) 
$ 2 = &Y u l + ^Y 2 -J, (99) 
12 

$ i = -y/2qv 1 ' + —q[v ij "C j + 6v ij 'v j - v j "C j + 6v p v tj ] (100) 
5 

-iV2\{i l v k " - 2 -^e J v l >"}e m , 

$ il = 2V2qv H ' + 2iV2qSymTF{C j v k " - \v ij "£ k )e jkl (101) 

-lSymTF{?>v i "£ l + 6v if v l + ^(6v ls v is ' + £'V S ")}. 
6 7 
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SymTF means take the symmetric trace-free part. 

When these results are substituted into the gravitational equation, 

[5T - ^T] = 2A^<& 1 ¥ 2 , (102) 

we see that the quadratic part of $1 drops out and plays no further role so that we 
could have written 

$i(T) = |. (103) 

Finally we have the solution of Eq. ( 11021) 

r 2 ° = -(TfT, (104) 
T = -T + T*Y° + T^., (105) 



where 



T e = -kj{-2V2if' + ~[{3£ 3 -£ j )lf j " + 3if j/ {7v j -v j ) (106) 
-3(f j - 3f V" + 3^'(5^ + v ej )]} 



T a = —kji^^^ -^SymTF[3(^ + + 3v i '(3v l -v l ) (107) 
+ '^[v is '(7v ls - v ls ) + (f" + tV'}}} 

+ t -^SymTF[(f ~ - \(? ~ V 

Before we go on to the evolution equation, Eq. (l69l) . we remark that these cal- 
culations were extremely long and rather ugly. They however were checked by doing 
them other ways; via the L version, Eqs.(l66l). ( 1671) and (1681) and then partially by the 
RTMaxwell versionP 

C. The Bondi Mass Aspect 

In this section we calculate the Bondi mass aspect, Eq. ffOl . 

^ = = + 2Lh(a)' + L 2 (a)- + 5 2 cr + a(a)\ (108) 

and then, from the I = (0, 1) harmonics, extract the energy-momentum four vector. 
Again the calculations are long and complicated with frequent uses of the Clebsch- 
Gordon expansions. 
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First, defining the harmonic coefficients by, 

^ = r + ry° + r^V4-, (109) 

and using 

r 2 ° = [-T + TYy + T^Y 2 %}T ' 3 , (110) 
we obtain, from Eqs. (l89aD . ( TT041) and (fTU5l) . that 

T = -2v^T [l + ^V + t;T + ^(eV J, + 2^V J ')] (111) 



T l = -6T v l + 2kq 2 v 1 ' + T ^^(v lj 'e +v j, £ lj +4v lj v j ) (112) 

5 

-kg 2 ^[(3f - e)v lj " + 3v lj '(7v j + v j ) - 3(f j - 3£ lj )v j " + 3v j, (5v lj + 7v lj )} 
-k q H[{t - C)v kl - - Wea, 



T ik = 12T v kl -^-v kl ' (113) 

3 

-r V2SymTF[v l '£ k + 2v l v k + -^{£ kl v ai + 2v kl v il )} 

+ ^^.SymTF{(f + tW" + v*(3i7* + 5v k ) + ^-[v is \7v ks + 5v ks ) 
6 7 

The remaining terms in the mass aspect, Eq. (11081) . i.e., 

2L5(a) ■ + L 2 (ay + 5 2 cf + a{a)' , 

are calculated from Eqs.fljl]), dZSI), (tHZD and (l89al) . 

Using this result and the expression for ip^ , we obtain 

^ = #(°) + ^rWyo + vMy&j, (114) 

with 

1 1 44 988 
#(0) = _2v / 2{T [l + -vV + (f V* + 2v'V)] + — C j v ij (115) 

-kq^'v' + ^-v ij 'v ij ]}, 
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tf« = -6{T v l - \kq 2 v 1 ' - To^A^'e + v j '^ + W) - ^?#^{ll6) 

+fcg 2 ^[(3f - f )#'" + 3u ,i/ (7v' + v j ) - 3(^ - 3£ lj )v>" 
+3^'(5# + 7v«)] + \ikq\f - C)v k ' - y (? j - ^)#T^}, 



^(ifc) = 12To ^ 1 - 1^!^' + 24^ fc (117) 
—T ^SymTF[v i '£ k + 2v l v k + — (f fc V + 2v kl v il )} 

+ ^^SymTF{(f + £ fc )lf" + #(3?J fc + 5i; fc ) 
6 

+y [^'(7^ + 5^ s ) + (f " + £ ks )v ls "} 



■(24) 2 v / 2. 



_^ mTF[ ^_j_^ ( 2^-e ) - i6iu w (e -rKn fe ], 

Finally, the condition 

# = (118) 
yields the reality conditions to be imposed on the harmonic components 

^(o)_^(o) = 0> ( lig j 

VI. The Bondi Energy-Momentum Four- Vector 

The relationship of the mass aspect, ^ rea i = ^ = -^p-, to the Bondi energy- 
momentum vector, 

P a = (Mc,P i ), (122) 

■„ I19I18I8I11IT2TT61 

*«ai = - ^PT« + ¥ l Y°, (123) 

By comparing with Eqs. (11231 ) with the final form of the mass aspect, Eq. (1114l ). - 
(which is now real) - the mass, M and momentum, P l can be found. 
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Using the relation 

e = e R +iCi, etc., 

leads to the kinematic relations between mass and momentum with the position, 
velocity and acceleration, quantities defined from the UCF, Eq. (l80l) . i.e., 



M = M {1 + ^(v^ r -v)v]) + ^(vUr-v^)) (124) 



he 2 

+^(fi4 + CM) - §^{44 + 4<4 + f(vM + vfvf)}, 



M v e R - g< - ™^*>[vg& - vfz + viei - vui + 4(44 - vim 



5G L H ^ H 1 15c 4 

"'■"'■'' U 3,-' L ^'< 5 



+3^4" + e&v? + 18^4 + 344] + Sei4 - ?e?4']W, 



and 



*<*> = ^4-^4' + 24^ (126) 

^ Mo WfK^ - + 2^4 - 2^14 + - 



^4 — c7 l ~ ris n ~ l si ■ — rt-rt — j-j 1 ij 

,kl„,il o„.kl„,il\ 



+2v R % - 2vfvf)} 

1 v 2(TV? %mTF[£*4 + 444 + 44 + ?(64'4 S - 4'4 S + &v*n 



(24) 2 v / 2 c __ r7ni7ir ,„,,,,k 



7c 

The reality conditions, Eqs. (1119D . (11201) and (11211) become 

= M [^4 + ^(*j& + ^ (127) 

+^(^4 - fi4) - ^Wi - 44 + f (4'4 - 4'4)L 
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5G J 15c 4 

and 



= — + -j^vf '- 24£* 



c° 

■^V^K^ + + 2^* + 2*44 + 2 -f(e R l vf + t# + 2v R l vf + 2v?v% 



+^{%vf -\$vfye jsk ] (129) 

(24) 2 V^2 771 TO.. IS ^ifcs ..ZS^fcs 



7c 

These 'kinematic' equations are of sufficient complexity that understanding of 
them is difficult or even problematic. It appears that perhaps the best way to proceed 
is to leave out for the time being some of the details - mainly the affect of the / = 2 
terms, e.g., v R , v 1 /) and their derivatives. 

We then obtain for the mass 

M = M {1 + ^(v R v R - v)v\) + - - 2 -^{vlv R + v«v}}, (130) 

the rest mass and kinetic energy term plus a series of unusual quadratic terms that 
do not appear to have any standard counter parts. If - as we argue later - we identify 

d l E = qC R and d l M = qd, 

as the electric and magnetic dipole moments and then arguing from the Kerr and 
charged Kerr solutions that we identify 

S' 1 ' = M cu\, 

as the spin dipole moment and its time derivative, we can interpret the extra terms 
in the mass as the interactions of these moments with themselves and each other and 
write: 

M = M o{1 + i(^<-^S^0 + i«4-^S'»S')} (131) 



9 

5 \ a E a E * U M U MJ- 
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.The linear momentum expression becomes 

P e = M v e R - + ^K^]'e lke , 

with the first term being the standard Mv, the second the contribution from radiation 
reaction and the last term is new. 

The first and second reality conditions become 

= M [v R v} + (vU) + vfC R )} - 2 -^Kv) - v}'v R ], (132) 

= M vf + ^vf - |![^<]'e ifce . (133) 

We see from Eq. (11330 that when the charge vanishes vf = and the spin S* is 
constant. From Eq. (11321) . S l is normal to v R . When q ^ 0, vf exponentially decreases 
and S* is asymptotically normal to v R . 

The third reality condition has the £f diverge exponentially at linear order 
though it is not clear what are the affects of the non-linearity. 

It appears that the reality conditions play the physical role of describing the 
evolution of the spin (or so-called 'magnetic moments'). 

We thus have found from the mass aspect and the reality conditions the following 
information: (a) there is a kinematic expression for the total mass, i.e., the bare 
mass dressed with the kinetic energy and spin-velocity coupled terms, (b) for the 
momentum there is the usual kinetic term with the classical radiation reaction term 
plus an electric/magnetic dipole coupling., (c) the first two reality conditions, though 
not fully understood, determine the behavior of the spin vector, (d) while the I = 2 
reality condition, also not well understood, determines the imaginary or spin part of 
the quadrupole term. 



VII. Dynamics: Equations of Motion 

The equations that we have been dealing with have essentially been relational, 
i.e., they just giving the relations between the different geometrical/physical quan- 
tities. Essentially all the variables have been determined in terms of the harmonic 
coefficients of the UCF and their derivatives, i.e., (v °, v \ v iJ ), etc. and the 

T . 

The final dynamical equation that determine the evolution of these quantities is 
the Bondi equation, ( 1691) . 

iff- = a "a + k(j)* 2 %°. (134) 

By the harmonic expansion of Eq. (11231 ). and using the expressions for 02° an d 
a' in Eq. ( 11341) . we obtain the mass and momentum loss equations, 
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1 Q9r 2 a 2 

P e ' = ^(^-4S J >^-^(^4-4V)% e . (136) 

8\/2^ / if jet if i*e/\ 32^ . fc?/ ii/ kit iif\ 
+ —^-( V R v R + V I V I ) + 15^4 ( U H U / _ V W fl 

The mass/energy loss equation fjl35j) has the following simple meaning: 

• If, in the first term on the right side, we make the physical identification of 
(v l R ,vf) with the derivatives of gravitational mass and spin quadrupole moments 
via 

(v R >,v?) = Y ^(QtLiQt^ (137) 
it is exactly the Bondi gravitational quadrupole energy loss. 

• For the second term, if we identify, as we mentioned earlier, (?(£#,£j) with the 
electric and magnetic dipole moments, it is precisely both electric and magnetic 
dipole energy loss. 

• If, in the third term, we make the identification of q(v R , v 1 /) with the derivatives 
of the electric and magnetic quadrupole moments via 

q(v%,v?) = ^^-{D% f ec ,D^' ag ), 

we have the electromagnetic quadrupole energy loss. 

The momentum loss equation , (11361) . also has a simple physical interpretation, 
namely it is the source recoil from the radiated gravitational and electromagnetic 
momentum. 

Our final task is to substitute our dynamic expressions for the mass and mo- 
mentum, Eqs. (11241) and (11251 into the mass and momentum loss equations, (11351) and 
(11361) . obtaining the equations of motion for what we define as the center of mass, 



M Q v% 



J e 



2g 2 e „ , 2g 2 



(138) 



3c 

{ l ^^[vfe R - vfi\ + 4$ - v}'& + A{vfv R - vfvi)) 

+ ^^[« + - 8 -^f[e R vT - 2&f + l2vfv R - 9vfvj 

+seiv j R " + eeM" + is« + *4'v?) + ^[^e/vf]'e ike }' 



3c 4 L 5 



384c 2 
5G 



v I J v^e ike + —^{v R v R +v J jV J j ) + 



fail. 



64q 2 
15c 7 



kj/ iji 
A V R V I 6ike 



26 



with the term J e depending on the I = 2 harmonics. We see a standard left side, 
mass times acceleration, equal to the classical radiation reaction term plus a coupling, 
(cross-products) between derivatives of the electric and magnetic dipole moments and 
the driving term coming from the different quadrupole moments coupled to themselves 
and the dipole moments. 

We point out, since we have been keeping terms only to second order, that we 
have left out a natural cubic term that comes from the mass loss equation. If we write 

p k = Mv k + 

we see that 

p k > = Mv k ' + M'v k + 

It is the term M'v k that we have omitted. From Eq. (11351) we have the cubic term 
that could have been included in the equations of motion: 

P» = M«»W k {^(iM+^^ 

(139) 

Actually this cubic term, without the gravitational contribution, (the first cubic term), 
is well known in classical electrodynamics^ 

Summarizing what has been found from the dynamical equations we see; (a) From 
the mass loss equations we have been able to identify or give the physical meaning 
to the £ R , £j, £ y , etc., in terms of different moments, (b) From the Kerr metric, we 
identified Mc^\ with the spin vector and from the charged Kerr metric (and also from 
the mass loss equation) we associate with the magnetic dipole moment. This leads 
immediately to the Dirac value of the gyromagnetic ratio, g = 2, (c) Finally we have 
the momentum loss equation and the equations of motion for the 'center of mass', 
£ R determining the motion of an gravitational/electromagnetic rocket. Without any 
attempt at model building we have automatically found the well-known radiation 
reaction force - classically leading to the unpleasant run-away behavior - but now 
with damping forces. Though it is not at all clear if the damping terms are sufficient 
to stop the run-away behavior, nevertheless if there was an existence proof for the 
exact type II solutions to asymptotically approach the Kerr-Newman metric that 
would be a strong indication for the suppression of the run-away affects. On the 
other hand it might well be that the type II metrics with a Maxwell field are in fact 
unstable. 

VIII. Conclusion 

In this work we have returned to the study of the important class of Algebraically 
Special Type II metrics and tried to give them physical content. Many attempts over 
the past forty or so years have been made to solve them - with very few exact solutions 
- most work being on approximations. Existence theorems for the pure vacuum case 
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(in the Robinson-Trautman subclass^) were given and much approximate work was 
done on them but there appears to have been little progress on the Maxwell version. 
There have been several attempts, 134 1 35 1 via approximations, to relate, at least, the 
pure RT metrics to the more general class of asymptotically flat Bondi metrics where 
one has an understanding of asymptotic energy and momentum. However, to our 
knowledge, little or nothing along these lines was done for the RTM equations.^ 

Here, with the help of the UCF, many of the ideas have been unified with a 
possible resolution of the physical content of these metrics. We have even given the 
relationship of these metrics to the Bondi metrics, via the UCF. 

Though later in this section we will speculate on some of the physical significance 
of our observations, we first summarize our general construction. 

In an arbitrary asymptotically flat space-times, we started with families of regu- 
lar asymptotic shear-free null geodesic congruences, described by its (stereographic) 
angle-field L(u,(,C) a t 3 + , that are generated by arbitrary complex curves in H- 
space. This function, L, and the associated curve are then completely determined by 
the requirement that the I = 1 harmonic of the asymptotic Weyl tensor component, 
if) , should vanish. This leads to a unique complex world-line and UCF, u — X(t, (, (), 
which then is expanded in spherical harmonics. This general procedure was then ap- 
plied to the analysis of the type II metrics - using severe approximations. When the 
harmonic coefficients of X(r, C, C) were used in the Bondi mass aspect (in particular, 
the Bondi energy-momentum four-vector) and the Bondi evolution equation, we were 
able to assigned to them 'natural appearing' physical meaning. 

• First, by separating the real and imaginary parts, we could identify the 1 — 0,1, 
harmonics and their derivatives, as the center of mass position and velocity 
vectors for the real parts and the specific spin-angular momentum for the imag- 
inary part. Then, from the mass aspect, we saw how the Bondi energy and 
momentum took on a kinematic meaning in terms of these kinematic variables. 

• From the Kerr metric and the charged Kerr metric, we identified Mc£} with 
the spin vector. From the mass loss equation, we associate £,)) with the 
electric and magnetic dipole moments. This leads immediately to the Dirac 
value of the gyromagnetic ratio, g = 2. Other multiple moments extracted from 
^( r )C;0) were identified via the energy loss theorem. 

• Finally, using the kinematic expression for the momentum, in the Bondi momen- 
tum loss equation we obtained equations of motion for the 'center of mass', de- 
termining the motion of what could be identified as a gravitational/electromagnetic 
rocket. 

Without any attempt at model building we automatically had the well-known 
radiation reaction force - which classically leads to the unpleasant run-away behavior 
but now it contains damping forces. Though it is not at all clear if the damping terms 
are sufficient to stop the run-away behavior, nevertheless if there was an existence 
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proof for the exact type II solutions to asymptotically approach the Kerr-Newman 
metric that would be a strong indication for the suppression of the run-away affects.^ 
Several questions concerning these results immediately come to mind. 

• In what space is this motion taking place? 

• In what sense could it be said that this motion is observable? 

• Are there any possible consequences that can be derived from these results? 

There are tentative answers to each. As we pointed out earlier, the motion can 
be thought of as taking place in either the space of the complex Poincare translation 
group or in H-space with the velocity norm coming from either the Minkowski or 
H-space metric. Though neither of these spaces are observable directly, in principle 
- but almost certainly not in practice - the motion could be deduced by observing 
the special unique class of asymptotically shear-free null rays, and then studying the 
I = 0, 1 coefficients of the harmonic decomposition of their angular distribution. In 
other words the complex curve is there (if these shear-free rays could be found) but 
in a strange form. 

The real physical meaning of these results certainly is not at all clear. They 
do contain basic well-known classical equations: Newton's 2nd Law, rest and kinetic 
energy, the usual expression for momentum, the electromagnetic dipole radiation for- 
mula, the Dirac value of the gyromagnetic ratio, the classical radiation reaction force, 
the gravitational quadrupole formula. In other words there is a strong suggestion that 
they do have meaning. If that were the case then they should be taken seriously and 
one should look for new physical results and predictions. For example, in Eq. (1131f ) we 
can see other contributions to the mass/energy (in addition to rest energy and kinetic 
energy) coming from the spin and from the electric charge. Though very small they 
are being predicted. 

Another interesting possibility arises from the radiation reaction force. This term 
in the equations of motion has been a source of considerable difficulty for classical elec- 
tromagnetic theoryp2E21 Almost all attempts to study the classical motion of charged 
particles interacting with their own field has led to the unphysical run-away behav- 
ior. The usual answer to this problem has been to invoke quantum or quantum field 
theory.^ Our suggestion or conjecture is that when our procedure is applied to any 
asymptotically flat Einstein-Maxwell space-time (with the two world-lines coinciding), 
assuming that the space-time is well behaved and approaches Reissner-Nordstrom or 
charged Kerr, then there will be a radiation reaction term in the center of mass equa- 
tions of motion whose influence will be suppressed by gravitational affects. There 
already are tentative results in this direction. This calculation is being planned for 
the near future. To do the calculation without the assumption of the coincidence of 
the two world-lines, though planned for the future, will be quite difficult. 

As a final comment we mention that though there has been much mathematical 
work done on CR structures and GR (see Sec. IV, Remark 3), it appears to us that 
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this is the first time that a CR structure has made relevant contact with physical 
issues. 
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